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Разработана и реализована в среде программирования Ое/р/і/ математическая модель частотно-регулируемого асинхронного 
двигателя. Учитывается влияние насыщения железа магнитопровода на параметры двигателя. Характер нагрузочного момента 
- активный или реактивный. Предусмотрена компенсация падения напряжения на активном сопротивлении обмотки статора. 
Учтена механическая характеристика тормозного устройства. Математическая модель позволяет получить характеристики асин- 
хронного двигателя при работе его от преобразователя частоты и может быть использована при проектировании частотно- регу- 
лируемого асинхронного двигателя. 

Во многих областях применения не регулируе- ческого процесса и могут иметь низкие энерге- 

мый привод уступает место, как более перспектив- тические показатели. 

ному, частотно-регулируемому асинхронному при- 2. Модернизированные асинхронные двигатели, 
воду [1]. В качестве частотно-регулируемых асин- созданные на базе серийных, но с учетом спе- 

хронных двигателей (ЧРАД) используют следую- цифики их работы от преобразователя частоты 

щие асинхронные двигатели: (ПЧ). При проектировании таких двигателей 

1. Серийные не регулируемые асинхронные дви- снимаются ограничения по кратностям пуско- 

гатели общего или специального назначения, вого тока, имеется возможность использовать 

рассчитанные на питание от сети с постоянной нестандартное номинальное напряжение. По 

частотой и напряжением. Такие двигатели не возможности, в модернизированных двигателях 

всегда удовлетворяют требованиям технологи- для обмотки ротора используют материал с бо- 
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лее высокой удельной электрической проводи- 
мостью, изменяют обмоточные данные обмот- 
ки статора и т.д. 

3. Узкоспециализированные асинхронные двига- 
тели, спроектированные с учетом законов упра- 
вления частотой и напряжением ПЧ, диапазо- 
ном регулирования частоты, характером на- 
грузки и других факторов, специфичных для 
конкретной области применения частотно-ре- 
гулируемого асинхронного привода. Стоимость 
таких двигателей существенно выше, чем се- 
рийных и модернизированных. 

При выборе существующих или при проектиро- 
вании модернизированных ЧРАД возникает 
необходимость в количественной оценке их харак- 
теристик при работе от ПЧ. Для решения такой за- 
дачи, особенно на стадии проектирования ЧРАД, 
целесообразно воспользоваться математическим 
моделированием частотно-регулируемого асин- 
хронного привода. 

Математическая модель частотно-регулируемо- 
го асинхронного привода должна удовлетворять 
следующим требованиям: 

1. Простота в реализации. 

2. Учет влияния насыщения железа магнитопро- 
вода на параметры асинхронного двигателя. 

3. Учет законов регулирования частоты и напря- 
жения ПЧ. 


4. Учет механических характеристик нагрузочного 
момента и тормозного устройства. 


Для математической модели ЧРАД использова- 
на математическая модель обобщенной машины в 
синхронной системе координат и,ѵ [2, 3]. Система 
координат и,ѵ позволяет наиболее просто изменять 
во времени, не зависимо друг от друга, амплитуду 
фазного напряжения и частоту ПЧ по любому за- 
данному закону: 


-ш ■, ■ Л 

сіш 

и =і -гн — + со ■ у/ ; 

хѵ хѵ х х т хи ’ 

°=Сѵ-с + ^+К 
, , . , . , сісо,. 

и -=м<«> +ТЧГ- 


( 1 ) 


где: и ш , и 5Ѵ - напряжения; і ш , і„, і п , і„ - токи, г„ г г - 
сопротивления; у/ !и , і//„, і// га , у/„ - потокосцепления 
обмоток статора и ротора обобщенной машины; 
М ш и М с ({) - электромагнитный и нагрузочный мо- 
менты; / Е - суммарный момент инерции, р - число 
пар полюсов двигателя; со х и со г - угловая частота 
вращения системы координат и,ѵ и ротора. 

Система координат и,ѵ вращается синхронно с 
магнитным полем обмотки статора с угловой частотой 


о) х =2я/ х , 


( 2 ) 


где / х - частота на выходе ПЧ. 


Систему дифференциальных уравнений (1) до- 
полним в матричной форме уравнениями потокос- 
цепления обмоток статора и ротора: 
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Полные индуктивные сопротивления фаз обмо- 
ток статора Ц и ротора Ь; 

— 'г. Ан > 

Іх = С т “Ь С , 

г а 2 п т ? 



где Ь Лп и Ь Лп - индуктивности рассеяния фаз обмо- 
ток статора и ротора с учетом насыщения, Ь я - вза- 
имная индуктивность. 

Для решения дифференциальных уравнений ( 1 ) 
воспользуемся численным интегрированием. С 
точки зрения устойчивости численного интегриро- 
вания целесообразно решение дифференциальных 
уравнений ( 1 ) относительно потокосцеплений об- 
моток статора и ротора. Введем обозначения: 

к /Ь ; к =1 /X ; к =1 /Ь , (6) 

гз г тз 1 зг 8 тз’ т т тз> у ' 


где Ь т =1;Ь,—Ѵ т - определитель матриц коэффи- 
циентов ур. (3, 4). 

При принятых обозначениях (6) из ур. (3, 4) вы- 
разим токи обмоток статора и ротора: 

Си — * к гз —Ц* ги ' к т , 

ку=У / „-К-У / п,-к т , 

/ 7 Г ( 7 ) 

С, =ѵ т -к зг -цг зи -к т , 
і — ш -к —ш ■ к . 

гѵ г гѵ хг т хѵ т ^ 

Подставим (7) в (1) и приведем математическую 
модель (1) к виду, удобному для численного инте- 
грирования: 

%.„,=[«„ +ч ■ѵ.у-к -(ѵ,, -к -Ѵп, -к )]■<*; 

= К - ч • - п ■ (ѵч • к - ѵч • к )] • , 0 . 

= Кч -ч)-ѵч -к -О/ч -к -к )]■<*; 

= Кч - ч ) • ѵч - к ■ (ѵч • к - ѵч ■ к )] • <*• 

Напряжения и ш и обмоток статора выразим 
через напряжение на выходе ПЧ 

«»=^х-со8 (%-%), 1 9) 

Чѵ=^,-5Іп(Ф„-Ф,). ] 

Амплитуда фазного напряжения ІІ Х является 
модулем изображающего вектора напряжений об- 
мотки статора. Амплитуда фазного напряжения 
ПЧ, фазовый угол изображающего вектора напря- 
жений ср и и угол поворота (р к системы координа- 
тных осей и,ѵ относительно оси а трехфазной об- 
мотки статора задаются начальными условиями. 
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При скалярном управлении ПЧ начальными 
условиями целесообразно задать ср= 0 и <р к = 0. В 
этом случае в момент времени 1= 0 ось и системы 
координат и,ѵ и изображающий вектор напряже- 
ния ІІ Х направлены по оси фазы а трехфазной об- 
мотки статора, в ур. (9): и,=Ѵ х , и, ѵ =0. Начальными 
условиями в момент времени /=0 задаются пото- 
косцепления ц/ ш , у/ т і//„, у/ п , обмоток, частота со г 
вращения ротора, частота/. ПЧ. 

При заданных начальных условиях на каждом 
шаге численного интегрирования рассчитываются: 
частота вращения со х (2); приращения потокосце - 
плений (8); потокосцепления обмоток 

+ ^ уг ю = У* + *'/'„* | 

Ѵп=У',, + с1у/ п ;, ) 


токи в обмотках статора и ротора (7); электромаг- 
нитный момент 

М = Щ- ■ ѵ ■ (ѵ/ -і —у/ • і ), 

эм о г \т ни лѵ т $ѵ 8и 


где пг= 3 - число фаз обмотки статора, р - число 
пар полюсов двигателя; приращение угловой ча- 
стоты вращения ротора 


г/й),. = р ■ 



■ сіі; 


частота вращения ротора и текущее время 
ю,. = й),. + г/©,. , 1 = 1 + сіі. 

По токам ур. (7) на каждом шаге численного 
интегрирования рассчитываются: модуль / изо- 
бражающего вектора токов обмотки статора 


фазовый угол (рі изображающего вектора токов об- 
мотки статора 

<Р, =(р„ -аг е№(і„/і„); 
коэффициент мощности 

со 5ср а =со5 (<р,.); 


активная мощность двигателя 

Р а= Т^лЛіС05ф а . 

При пусках ЧРАД с предельным нагрузочным 
моментом с уменьшением времени регулирования 
і КІ и, особенно, с увеличением напряжения ком- 
пенсации II тч (рис. 2) токи в обмотках могут превы- 
сить номинальные в два и более раза. ЧРАД может 
выполняться с закрытыми пазами ротора. Для та- 
ких двигателей при токах, начиная с номинального 
и выше, из-за насыщения потоками рассеяния мо- 
стика закрытых пазов ротора, пазовое рассеяние 
обмотки ротора уменьшается. Короткозамкнутая 
обмотка ЧРАД может выполняться из меди 14]. 
Предельная верхняя частота может составлять 
100 Гц и выше при номинальной частоте/=50 Гц. В 
связи с этим, в математической модели целесооб- 
разно предусмотреть влияние насыщения железа 
магнитопровода потоками рассеяния на индуктив- 


ности рассеяния обмоток статора Ь л и ротора Ь аЪ а 
также влияние вытеснения тока в стержнях ротора 
на активное сопротивление г г обмотки ротора. Это 
позволит моделировать «жесткие» пуски, при кото- 
рых учет влияния насыщения и вытеснения тока 
наиболее актуален. 

Индуктивности рассеяния с учетом насыщения 
предлагается рассчитывать в функции тока обмот- 
ки статора, а активное сопротивление ротора с уче- 
том вытеснения тока - в функции частоты тока об- 
мотки ротора. На основании данных проектирова- 
ния двигателя считаются известными: не насы- 
щенные Ь л , и насыщенные Ц р , Ь 2р при скольже- 
нии 5=1 индуктивности рассеяния обмоток статора 
и ротора; пусковой ток І рХ обмотки статора и актив- 
ное сопротивление г 2р обмотки ротора при сколь- 
жении 5=1. Для учета влияния насыщения опреде- 
ляется разность между не насыщенными и насы- 
щенными значениями индуктивностей рассеяния 

^1 = Аті _ Ар’ ^А = Аг2 _ А ; р- 

В функции тока / рассчитываются насыщен- 
ные значения индуктивностей рассеяния обмоток 
статора и ротора 



где І кр - критический ток обмотки статора, ниже 
которого насыщение магнитной цепи не оказывает 
заметного влияния на индуктивности рассеяния 
обмоток. Величина этого тока определяется из 
предварительных расчетов механической характе- 
ристики двигателя. В случае закрытых пазов на ро- 
торе целесообразно принять ток І кр равным номи- 
нальному току обмотки статора. Для учета насыще- 
ния достаточно на каждом шаге интегрирования 
выполнить 3-4 последовательных приближений. 

Для учета влияния вытеснения тока в стержнях 
обмотки ротора на каждом шаге интегрирования 
рассчитываются скольжение 
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и коэффициент изменения активного сопротивле- 
ния ротора 


к = 0,2 + 


2 Р 


У/? 


- 0,2 


Ѵ аЬ Чл), 


где/ - базовая частота, на которую рассчитан дви- 
гатель, й),=2я/ - угловая частота, г 2 - сопротивле- 
ние обмотки ротора без учета вытеснения тока. 

Если к п > 1, то активное сопротивление обмотки 
ротора с учетом вытеснения тока рассчитывается 
по формуле г=к г ;г ъ иначе г=г г . 
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Математическая модель ЧРАД реализована в 
среде программирования БеІрЫ. Исходными дан- 
ными к моделированию являются номинальные 
данные и параметры двигателя. 

Номинальные данные: номинальная мощность 
двигателя Р„, номинальное фазное напряжение ІІ нф , 
частота сети/, число пар полюсов р, коэффициент 
полезного действия г) н , коэффициент мощности 
со 8<р, номинальное скольжение кратность пу- 
скового тока /, кратности пускового и критиче- 
ского моментов М р1 , М к1 , суммарный момент инер- 
ции ротора и приведенного к валу двигателя меха- 
низма /. 

Параметры двигателя учитывают индуктивные и 
активные сопротивления, соответствующие номи- 
нальному (без учета насыщения железа и вытесне- 
ния тока) и пусковому режимам работы двигателя (с 
учетом насыщения железа и вытеснения тока при 
скольжении 5=1). Параметры номинального режима 
работы: индуктивное сопротивление рассеяния об- 
моток статора х л и ротора х а2 , активное сопротивле- 
ние обмотки статора г. и ротора г„ индуктивное со- 
противление взаимной индукции х т коэффициент 
насыщения магнитной цепи в режиме холостого хо- 
да к ш . Параметры пускового режима: индуктивные 
сопротивления рассеяния х Лр и х Лр - обмотки стато- 
ра и ротора, активные сопротивление обмотки рото- 
ра г 1р . 

В программе математического моделирования 
нагрузочный момент М с (і) считается постоянным, 
но может быть задан в зависимости от частоты вра- 
щения ротора со,. активным, рис. 1, а, или реактив- 
ным, рис. 1, б. 



МДі) 
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Рис. 1. Механическая характеристика нагрузочного момента 


Величина нагрузочного момента задается в от- 
носительных единицах коэффициентом 
к ш =М с (і)/М н , где М н - номинальный момент двига- 
теля. В случае реактивного момента коэффициент 
к ж принимает знак направления частоты вращения 
ротора, т.е. к ж =к тс $щп(со г ). 

В программе математического моделирования 
предусматривается возможность регулирования 
частоты/ и напряжения ІІ Х ПЧ по линейному зако- 
ну, рис. 2, с компенсацией напряжения. 

Напряжение и частота ПЧ представлены на 
рис. 2 в относительных единицах (о.е). За базовые 
значения напряжения и частоты приняты устано- 
вившиеся значения напряжения и частоты. 

В исходных данных задаются времена: - регу- 
лирования по линейному закону частоты и напряже- 


ния ПЧ при разгоне двигателя; і т - работы двигате- 
ля при установившихся значениях частоты /, и на- 
пряжения ІІ уя ПЧ; 4, - останова двигателя, за кото- 
рое частота и напряжение ПЧ изменяются по линей- 
ному закону до своих критических значений / г и \] ь . 



Рис. 2. Диаграммы напряжения и частоты ПЧ: 1) напряжение 
4 , 2) частота к 


При активном характере нагрузочного момента 
задается 4 - время, в течение которого двигатель 
находится в заторможенном состоянии и после ко- 
торого происходит срабатывание тормозного 
устройства с последующим растормаживанием 
двигателя. При критических значениях частоты/,, и 
напряжения ІІ кг происходит срабатывание тормоз- 
ного устройства с последующим торможением дви- 
гателя и отключением его от ПЧ. 

Программа математического моделирования 
предусматривает возможность задания начальных 
(/о, / (при 1=0) и установившихся / й значений 
напряжения и частоты ПЧ. За время регулирова- 
ния 4 г частота и напряжение ПЧ изменяются от на- 
чальных до установившихся значений по линейно- 
му закону. 

Имеется возможность задания только началь- 
ной и установившейся частоты ПЧ. В этом случае за 
время і щ частота изменяется по линейному закону 

І гвё 


При выполнении условия (? ге8 +4„)>/>4 г прини- 
мается/^,. За время 1>(4 г +4ѵ) частота изменяется 
также по линейному закону 


Г 

Ух У узі ^ V ге§ гаѵ /’ 


Напряжение ПЧ является функцией частоты 
11= Д/) и с учетом напряжения П шм компенсации 
вычисляется по формуле 


ѵ =Ѵ2- 


II... 


V .-V 

нф КС 


/і 


■У 


где ІІ„ф и / - номинальные значения фазного на- 
пряжения и частоты напряжения двигателя. 

Тормозное устройство включается (рис. 2) при 
разгоне двигателя в момент времени 4, при тормо- 
жении - в момент времени ?>/. В исходных данных 
для математической модели тормозного устройства 
задаются: 4 - время растормаживания и М тм — мак- 
симальный тормозной момент. Тормозной момент 
является функцией времени и частоты вращения 
ротора. В функции времени тормозной момент из- 
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меняется за время ( !Г по линейному закону - от мак- 
симального значения до нуля в режиме расторма- 
живания и от нуля до максимального значения в ре- 
жиме торможения. В функции частоты вращения 
ротора тормозной момент носит реактивный харак- 
тер, принимая знак направления вращения ротора. 

Апробация изложенной выше модели была прове- 
дена нами на базе кранового восьмиполюсного ЧРАД 
4МТКД280Ь8 с короткозамкнутым ротором мощно- 
стью 90 кВт [4]. Двигатель спроектирован и изготовлен 
ОАО «СКВ Сибэлектромотор», г. Томск, на базе крано- 
вого двигателя серии 4МТМ280Ь8 с фазным ротором 
и предназначен для частотно-регулируемого асин- 
хронного привода механизма подъема башенного кра- 
на. В качестве примера на рис. 3 представлены резуль- 
таты моделирования. Подъем груза осуществляется в 
двигательном режиме, спуск груза - в генераторном. 

За базисные момент и ток приняты номиналь- 
ные значения момента и тока двигателя. За базис- 


ные напряжение и частоту приняты установивши- 
еся значения частоты и напряжения ПЧ. За базис- 
ную частоту вращения принята синхронная часто- 
та вращения ротора, соответствующая установив- 
шейся частоте ПЧ. 

Номинальные данные и параметры двигателя 
4МТКД280ІД являющиеся исходными данными к 
моделированию (рис. 3), представлены в [5]. Со- 
гласно техническому заданию (ТЗ) двигатель дол- 
жен обеспечивать постоянство момента в диапазо- 
не частот от 5 до 50 Гц и постоянство мощности на 
валу двигателя при частотах от 50 до 100 Гц. По ТЗ 
предельный момент нагрузки М=\,2М„. Результа- 
ты моделирования (рис. 3) соответствуют предель- 
ному моменту нагрузки, установившейся частоте 
/ м =5 Гц, напряжению компенсации (7 гал ,=10 В, вре- 
мени регулирования напряжения и частоты при 
пуске двигателя і к = 4 с, механической характери- 
стике тормозного устройства - М т =\,1М„, 1= 0,4 с. 




Рис. 3 . Временные зависимости (во.е.): 1) фазный ток статора, 2) электромагнитный момент, 3) коэффициент мощности, 4) ча- 
стота и 5) фазное напряжение на выходе ПЧ, 6) частота вращения ротора 
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На основании математического моделирования 
установлено следующее. 

1. Принятый закон регулирования частоты и на- 
пряжения ПЧ (при напряжении компенсации 
[І кш = 1 0 В) обеспечивают надежный подъем и 
спуск груза двигателем 4МТКД280Ь8 в соответ- 
ствии ТЗ. 
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